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SUMMARY 

Probabilistic composite micromechanics methods are developed that simu- 
late uncertainties in unidirectional fiber composite strengths. These methods 
are In the form of computational procedures using composite mechanics with 
Monte Carlo simulation. The variables for which uncertainties are accounted 
Include constituent strengths and their respective scatter. A graphite/epoxy 
unidirectional composite (ply) is studied to Illustrate the procedure and its 
effectiveness to formally estimate the probable scatter in the composite uniax- 
ial strengths. The results show that ply longitudinal tensile and compressive, 
transverse compressive and Intralaminar shear strengths are not sensitive to 
single fiber anomalies (breaks, interfacial disbonds, matrix microcracks); how- 
ever, the ply transverse tensile strength is. 


INTRODUCTION 

The analysis of composite structures requires reliable predictive models 
for material properties and strengths. However, the prediction efforts have 
been complicated by Inherent scatter in experimental data. Since uncertainties 
in the constituent properties, fabrication variables, and Internal geometry 
would lead to uncertainties in the measured composite properties, the question 
arises: 

How much of the "statistical" scatter of experimentally observed composite 
properties can be explained by reasonable statistical distribution of 
input parameters (primitive variables) in composite micromechanics and 
laminate theory predictive models? 

In order to answer this question, a study was conducted to develop a com- 
putational simulation procedure for probabilistic composite micromechanics 
(ref. 1). Application of this approach for uniaxial thermal and mechanical 
properties is summarized in reference 2. The objective of the present paper 
is to describe this type of micromechanics for fiber composite uniaxial 
strengths and present typical results obtained therefrom. The computational 
simulation is performed using ply substructuring with an existing computer 
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code (ref. 1) for composite mechanics and In conjunction with Monte Carlo simu- 
lation. The scatter In the constituent strengths Is selected from anticipated 
respective probabilistic distributions. 


COMPUTATIONAL SIMULATION FOR PROBABILISTIC COMPOSITE MICROMECHANICS 

In this section, the formal approach to computationally simulate probabi- 
listic composite micromechanics Is summarized. 


De termini stlc/Probabl 1 1st 1c Model 

The model commonly used for deterministic composite mechanics Is based on 
the calculation of properties of the basic unit of an orthotropic ply. The 
layup geometry Is then used In laminate equations to calculate composite prop- 
erties as shown schematically In figure 1(a). In the probabilistic simulation 
however, the basic unit Is taken as a subply (ply substructuring) which con- 
sists of only a single fiber-matrix. Deterministic micromechanics theory 
(refs. 3 and 4) Is used to predict the properties of the assumed orthotropic 
subply. The probabilistic aspect Is Introduced by representing the scatter In 
the fabrication variables and constituent material properties. Probable fiber 
misalignment within the ply are then used In the laminate theory equations to 
predict ply properties. This substructuring of the ply represents a novel 
attempt at characterization of fiber composite material properties based on 
probabilistically distributed constituent properties, Individual fiber misal- 
ignment and fabrication process (primitive) variables as shown schematically 
In figure Kb) . 

Ply substructuring in conjunction with composite mechanics Is particularly 
well-suited to the probabilistic description of fiber composite material prop- 
erties. The micromechanics and laminate theory equations can be used to calcu- 
late ply properties at any number of points In a ply. This approach provides 
a rational procedure for composite material property assessment because It 
evaluates ply behavior as the result of a series of random events (uncertain- 
ties In the primitive variables) which occur at the Intraply or micromechanics 
level. 


Composite Mechanics 

The probabilistic simulation Is performed by considering the ply as an 
assembly (equivalent laminate) of 15 subplles. The composite mechanics used 
In the simulation Is that available In the Integrated Composite Analyzer 
(ICAN) (ref. 5), which Is a computer program for comprehensive linear analysis 
of multilevel fiber composite structures (ref. 5). The program contains the 
essential features required to effectively design structural components made 
from fiber composites. It now represents the culmination of research conducted 
since the early 1970's at the National Aeronautics and Space Administration 
(NASA) Lewis Research Center (LeRC), to develop and code reliable composite 
mechanics theories. This user friendly, publicly available code Is depicted 
schematically In figure 2 and Is described In detail In reference 5. 
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Probabilistic Simulation - Monte Carlo Methods 

Complicated probabilistic events can be simulated by a variety of methods 
generally referred to as Monte Carlo methods (ref. 1). The term refers to 
that branch of mathematics concerned with numerical experiments on random num- 
bers. Since the advent of high speed computers, they have found extensive use 
In most fields of science and engineering, in analyzing many physical processes 
of a probabilistic nature, or where physical experimentation is limited or not 
feasible. 

A Monte Carlo simulation refers to the procedure where a single computa- 
tional simulation is performed by randomly assigning a value to an Independent 
random variable in a chosen model, and observing the dependent variable at the 
conclusion of the process being modeled. A Monte Carlo simulation is composed 
of n such Independent simulations. When n is sufficiently large, the 
observations will yield a statistically meaningful description of the physical 
problem. 

The form of Monte Carlo simulation used in the present investigation is 
as follows: 

(1) Define the model by assuming that: 

(a) It represents the composite mechanics 

(b) It is formulated In terms of primitive variables 

(c) It has probability distributions for the scatter in each primitive 
variable 

(2) Use the computer and random sampling techniques to select values of 
the primitive variables from their respective distributions. 

(3) Calculate dependent response variables using the model. 

(4) Replicate the experiment, each time with a new set of randomly sampled 
input values. 

(5) Use appropriate statistical methods to calculate probability 
distribution of the properties of Interest. 

(6) Estimate regression parameters for the assumed model. 


Computational Simulation Procedure 

To perform the computational simulation, a computer code was developed to 
couple ICAN and an available statistical analysis code (ref. 1). The logic 
diagram for this code is shown in figure 3. The steps are as follows: 

(1) Select values for the primitive variables for each subply from their 
respective assumed probabilistic distributions (fig. 4): 

(a) Normal - constituent elastic properties and fiber volume ratio, 
and fiber mi sal ignment 
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(b) Welbull - constituent strengths 

(c) Gamma - for void volume ratio 

Fifteen different sets (one for each subply) are generated where the means and 
scatter ranges for the primitive variable were those typical for AS graphite/ 
epoxy composite (tables I and II). 

(2) Enter these values as Inputs into ICAN. 

(3) Run ICAN and retrieve and store ICAN output for desired ply 
properties . 

(4) Repeat the process n-times where n is sufficiently large 

(50, herein) to provide data repeatability with an acceptable level 
of confidence. 

(5) Process the stored output using statistical analysis for cumulative 
probability distributions, confidence levels and significance. 


RESULTS AND DISCUSSION 

Results obtained by using the computational procedure described previously 
are presented and discussed in this section. The results are for ply uniaxial 
strengths, longitudinal tension (Sanj), longitudinal compression (S^iiq). 
transverse tension (S^22 T ^ * transverse compression (S&22C* ancJ Intralaminar 
shear (Sq,i 2S ) - ^e results are presented in graphical form where the ply uni- 
axial strength is plotted versus fiber volume ratio for ranges of scatter as 
represented by the Wei bull shape parameter (a). 


Longitudinal Tensile Strength (Sq.ht> 

The Influence of the scatter in fiber tensile strength on the ply longitu- 
dinal tensile strength Is shown in figure 5. Recall that the greater the value 
of a, the smaller the scatter. It Is seen that the spread in in- 

creases with Increasing fiber volume ratio. If we assume that the fiber 
strength scatter Increases with increasing fiber volume ratio due to fiber con- 
tact (abrasion) damage during processing, the following is deduced from the 
figure: the scatter in the ply longitudinal strength will also Increase. For 

example, we can see approximately from the figure, the anticipated range for 
SftllT to be 132 < Sj^nj/Sfy < 144 ksi for a mean fiber volume ratio of 0.5 
and a mean fiber strength (Sfy) of 400 ksi (table I). This range is lower 
than SjuiT = 200 ksi which is estimated using deterministic composite microme- 
chanics (table III). 

It is Important to recall that the probabilistic ply strength was pre- 
dicted by assuming that the ply fractures when the weakest fiber (subply) 
through its thickness fractures. Since the deterministic composite microme- 
chanics value is considerably higher (about two times), the conclusions are 
that: (1) a single fiber break is not sufficient to fracture the ply, (2) fi- 

ber load redistribution must take place, and (3) several fibers through the 
ply thickness must break prior to ply fracture. This is consistent with what 
is common knowledge in the composites community. It indirectly demonstrates 
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the substantial fracture toughness Inherent in composite longitudinal tensile 
strength relative to isolated fiber breaks. These Isolated fiber breaks have 
negligible influence on ply strength/fracture. As a side comment, the probable 
number of fiber breaks prior to ply fracture can be computationally simulated 
by accounting for fiber progressive fracture with simultaneous load redistribu- 
tion. The computational procedure will be analogous to that in CODSTRAN (Com- 
posite Durability Structural Analysis) for deterministic progressive composite 
fracture (ref. 6) but applied to single fiber breaks. 


Longitudinal Compressive Strength (Sjmc) 

The influence of the scatter in fiber compressive strength on ply S^hq 
strength is shown in figure 6 for the same a ranges as for The 

Important observations from this figure are (1) the curves peak at about 
0.5 FVR and then decrease, and (2) the curves for each a do not remain order- 
consistent but cross over. Both of these occur because of the four different 
fracture modes that are assumed to Induce ply longitudinal compressive frac- 
ture (ref. 4). It is also observed that the spread in S^ic changes for the 
different a's with increasing FVR. For example, the scatter between a 
equals 20 and 10 is greater for FVR 0.4 than for 0.6 and It is outside this 
range. One Important conclusion from these observations is that it would be 
very difficult to make consistent conclusions from experimental data for 

s anc- 

The range for S allc at 0.5 FVR from figure 6 is estimated between 117 
and 122 ksi which Is lower compared to the micromechanics estimate of about 
165 ksl. Again, the explanation is that a single fiber fracture does not 
cause ply fracture. And the ply Is relatively insensitive to single fiber 
anomalies. This Is a significant finding in view of the prevailing contention 
that single fiber anomalies are detrimental to ply longitudinal compressive 
strength. 


Transverse Tensile Strength (Sjj,22T> 

The influence of the scatter of the matrix tensile strength on SQ-22T Is 
shown in figure 7 for the same values in a as for Sjujy and Sgjic- The 
observations to be noted are: (1) the curves for the different a's remain 

order-consistent; that is, no cross-over occurs, (2) the curves indicate that 
S&22T continuously decreases with increasing FVR, and (3) the spread in the 
scatter decreases as the FVR Increases. 

The range in Sq, 22T between 9.0 and 10.4 ksl, which is lower than 
12 ksl (table III) predicted by deterministic composite micromechanics. How- 
ever, it is within the range of experimentally observed data for this strength 
(ref. 4). The authors consider this an important finding because it suggests 
that transverse tensile strength may be strongly Influenced by single fiber 
disbonds or microcracks. 


Transverse Compressive Strength ( Sg,22C^ 

The influence of the scatter of the matrix compressive strength on S& 22 C 
is shown in figure 8. The Influence is similar to that for S&22T- The 
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observations made for S^22T are applicable to Sg.22C as wel1 except for the 
specific values comparisons. The scatter In Sq, 22C is between 20.6 and 
23.8 ks 1 at 0.5 FVR which is considerably smaller than 27 ksl (table III) pre- 
dicted by deterministic composite micromechanics which Is also less than the 
range of experimental data of about 30 to 35 ks 1 (ref. 4), 

The authors attribute this difference to indicate that the ply transverse 
compressive strength is not sensitive to single fiber anomalies as is the ply 
transverse tensile strength. 


Intralaminar Shear Strength (Sq,]2S^ 

The influence of the scatter of the matrix shear strength on Sq,]]s is 
shown in figure 9. The Influence is similar to those for S^22T and Sq, 22C 
as would be expected since all of these are matrix controlled properties. 

There is some difference in the variation of scatter with FVR. The greatest 
spread is at about 0.5 FVR for Sqj 2S> whereas it progressively decreases with 
increasing FVR for the other two. 

The range in scatter for S^i 2S at 0-5 FVR is between 11.8 and 14.6 ksl 
compared to 10 ksi estimated by deterministic composite micromechanics. The 
authors Interpret this good comparison to indicate the following two signifi- 
cant points (1) the In situ matrix shear strength Is not influenced by the fab- 
rication process and (2) the ply Intralaminar shear strength Is not sensitive 
to single fiber anomalies. 


GENERAL DISCUSSION 

The scatter in the ply uniaxial strengths described and discussed was lim- 
ited to that influenced by respective scatter in fiber and matrix strengths. 
Other factors Influence the scatter as described in references 1 and 2. These 
factors Include scatter caused by single fiber anomalies In strength, interfa- 
clal disbonds and matrix microcracks. 

The authors consider the probabilistic simulation described in this paper 
as an Illustration of what can be done to quantify the uncertainties associated 
with the numerous factors that Influence ply uniaxial strengths. They do not 
consider the simulation to be complete nor the graphical results presented and 
the respective numerical values discussed as absolute. The authors strongly 
believe, however, that this is a rational approach to formally represent uncer- 
tainties associated with various factors that influence ply uniaxial strengths. 
The results obtained to date demonstrate the authors' contention. 

Regression results presented in reference 1, but not summarized here due 
to space limitation, indicate that factors influencing uniaxial different ply 
strengths are important in specific ranges of FVR and that no generalizations 
can be made at this time. 

The authors hope that the description, results and discussion summarized 
herein, will stimulate other investigators to pursue probabilistic representa- 
tion of composite uniaxial strength behavior beyond the longitudinal tensile 
strength which has extensively been Investigated over the years. 
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SUMMARY OF RESULTS 


The Important results of an Investigation to computationally simulate the 
probable scatter in composite uniaxial strengths as Influenced by scatter in 
respective fiber and matrix strengths are summarized below. 

1. A computational procedure has been described for probabilistic compos- 
ite micromechanics for uniaxial strengths. 

2. The scatter range In the uniaxial strengths Is represented In terms of 
the Welbull shape function In the respective constituent material strengths 
(fiber and matrix) for different fiber volume ratios. 

3. Comparisons with respective deterministic mean values and corresponding 
experimental data Indicate that ply longitudinal tensile, longitudinal compres- 
sive, transverse compressive and intralaminar shear strengths are not sensitive 
to single fiber anomalies. However, the ply transverse tensile strength Is. 
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TABLE I. - COMPUTATIONAL PROCEDURE INPUT 


DATA-FIBER PROPERTIES 


Probab 

i 1 i sti c 

Case 2 

Case 3 



0.0 


5.0 

10,0 

0.5 

0,5 

0.1 


0.3 

0.3 

0.03 

0.05 

con 

31 000 

HU 

3 000 

2 000 

2 000 

100 

200 

2 000 

2 000 

100 

200 


1 000 


100 


0.2 

HD 

0 

0.2 

0.2 

0 

0 

5.6 

5.6 

0 

0 

400 

400 

20 

10 

400 

400 

20 

10 


Input 


Oetenni ni st i c 


Case 1 


Theta, deg 

P 

a 

FVR 

P 

a 

WR 

X 

K 

E fll* ksi 
P 
a 

E f22 . ksi 

n 

a 

Gf 12* ^ 

P 

a 

G f23 , ksi 
P 
a 

v f 12 
P 
a 

a f11> Ppm/°F 
P 
a 

22 * PP m/ °F 

P 

a 

S f t » ksi 


S fC , ksi 


0.0 


0.50 
0.01 
31 000 
2 000 
2 000 
1 000 
0.2 
0.2 
5.6 
400 
400 
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TABLE II. - COMPUTATIONAL PROCEDURE INPUT DATA 


MATRIX PROPERTIES 


Input 

Determi ni sti c 

Probabi 1 i sti c 

Case 1 



E m , ksi 

500 


lIMPBIlilM 

P 

— 

500 

500 

a 

— 

25 

50 

v m 


0.35 


M 

— 

0.35 


a 

— 

0 


“m> Ppm/°F 

36 

— 

HI 

P 

— 

36 

36 

a 

— 

0 

0 

s mT' ks1 

15 



— 

P 

— 

15 

15 

a 

— 

20 

10 

S mC , ksi 

35 

— 

— 

is 

— 

35 

35 

a 

— - 

20 

10 

s mS ’ ksi 

13 

— 

— 

0 

— 

13 

13 

a 

— 

20 

10 


TABLE III. - DETERMINISTIC 
PLY PROPERTIES 
[0.5 Fiber volume ratio.] 


























ORTHOTROPIC PLY 
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(a) CONVENTIONAL — MULTI-FIBER PER PLY. 
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SUBPLY 



(b) PLY SUBSTRUCTURING - SINGLE FIBER PER PLY. 
FIGURE 1. - PLY SUBSTRUCTURING ANALOGOUS TO LAMINATE DE- 
COMPOSITION, SUBSTRUCTURING. 


COMPONENT 


/ 


/ 


/ 


LAMINATE 


\ 


\ 


\ 


\ 


\ 


\ 



\ 


\ 


\ 


LAMINATE 


COMPOSITE 

MICROMECHANICS 

THEORY 


PLY 


COMPOSITE 

MICROMECHANICS 

THEORY 


/ 


/ 


/ 


/ 


UPWARD 

INTEGRATED 

OR 

"SYNTHESIS" 


\ 


CONSTITUENTS 
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FIGURE 2. - ICAN: INTEGRATED COMPOSITIES ANALYZER. 




MONTE CARLO SIMULATION OF UNCERTAINTIES IN PRIMITIVE VARIABLES 


QUANTIFICATION OF PLY PROPERTY UNCERTAINTIES 


CALL 

SUBROUTINE 

COPY 


WRITE 
INPUT TO 
NEW 
FILE 


CALL 

SUBROUTINE 

ICANMN 




FIGURE 3. - PROBABILISTIC COMPOSITE MICROMECHANICS COMPUTATIONAL PROCEDURE LOGIC DIAGRAMS. 
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CUMULATIVE FREQUENCY 



RANGE (E -01) RANGE (E -01) RANGE (E 00), Ksr 

HISTOGRAM HISTOGRAM HISTOGRAM 



RANGE (C -01) RANGE (E 01) RANGE (E 00), ksi 

CUMULATIVE DISTRIBUTION CUMULATIVE DISTRIBUTION CUMULATIVE DISTRIBUTION 


NORMAL 


ga™ 


WE I BULL 


FIGURE - ASSUMED PROBABILITY DENSITY FUNCTIONS FOR UNCERTAINTIES IN PRIMITIVE VARIABLES. 



FIBER VOLUME RATIO 


FIGURE 5. - FIBER TENSILE STRENGTH SCATTER EFFECTS ON LONGITU- 
DINAL TENSILE STRENGTH. 


12 



FIBER VOLUME RATIO 


FIGURE G. - FIBER COMPRESSIVE STRENGTH SCATTER EFFECTS ON LONG- 
ITUDINAL COMPRESSIVE STRENGTH. 



FIGURE 7. MATRIX TENSILE STRENGTH SCATTER EFFECTS ON TRANSVERSE 
TENSILE STRENGTH. 
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A o = 20 
□ a = 15 
V a = 10 
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